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Abstract 

We propose a novel mechanism for the origin of non-Gaussian tails 
in the probability distribution functions (PDFs) of local variables in 
nonlinear, diffusive, dynamical systems including passive scalars ad- 
vected by chaotic velocity fields. Intermittent fluctuations on appro- 
priate time scales in the amplitude of the (chaotic) noise can lead to 
exponential tails. We provide numerical evidence for such behavior in 
deterministic, discrete-time passive scalar models. Different possibili- 
ties for PDFs are also outlined. 

PACS numbers: 02.50Ey,05.40+j,47.27Qb 
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The explanation of the form of the probability distribution of the fluctu- 
ations in local variables in extended, dissipative, dynamical systems poses an 
interesting challenge. In various turbulence experiments ||, ||, |J, Probabil- 
ity Distribution Functions (PDFs) of velocity gradients and passive scalars 
(temperature) are observed to be non-Gaussian; in particular, exponential 
tails have been seen whereas Gaussian distributions might be expected if one 
naively invokes the central limit theorem. Theoretically, a variety of expla- 
nations including an ingenious phenomenological model based on a nonlinear 
Fokker-Planck equation || and methods that rely on obtaining closure and 
moment balance relations |7| have been proposed to explain exponential 
tails in the passive scalar problem. The moment balance approach has been 
used to provide good fits to experimental data . In addition, several eddy 
diffusive models that yield exponential tails have also been studied 0. 

In this letter we explore a new mechanism in which the nature of the tem- 
poral correlations of the fluctuations that couple to a diffusive field plays a 
crucial role. We focus on passive scalar models and show that non-Gaussian 
PDFs including exponential and stretched exponential behaviors can arise 
when local variables are coupled to fluctuations whose amplitude varies ran- 
domly on a time scale comparable to the intrinsic diffusive time scale. In 
experimental fluid systems, exponential tails could, therefore, arise if the 
velocity fluctuations exhibit intermittent behavior in their magnitude on rel- 
evant time scales. 

We consider the diffusion of the temperature field 0(r,t) advected by a 
flow characterized by the velocity V in the presence of an externally-imposed 
mean temperature gradient (3. This is described by the passive scalar equa- 
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tion 

^ = uV 2 9 -V-V9. (1) 
dt v 1 

The coefficient v is the effective thermal diffusivity and the velocity is as- 
sumed to be incompressible: V • V = 0. 

The behavior of (1) in the presence of prescribed random velocity fields 
can be studied numerically directly in the continuum. However, the case of 
a velocity field arising from a deterministic chaotic dynamics is computa- 
tionally prohibitive if one uses Navier-Stokes equations. Instead we consider 
2-dimensional coupled map lattice models in discrete time and space which 
can be viewed as coarse-grained versions of the continuum equation with the 
lattice spacing of the order of the correlation length of the velocity fluctu- 
ations. The scalar 8(i) is defined on each site i of a L x L square lattice 
with a mean gradient (3 along the y direction and periodic boundary condi- 
tions in the x direction. The system evolves synchronously in discrete time 
n according to 

e n +i(i) = e n (i) + vV 2 L e n (i) - v n (i) ■ v L e n (i) (2) 

where Vl is a symmetric, lattice gradient. The incompressibility of the two- 
dimensional velocity field V n {i) = (u n (i), v n (i)) is enforced by obtaining it 
from a stream function ip n {i): u n {i) = •yVyip n (i);v n (i) = — •y'V x tp n (i) . The 
velocity vanishes at both y = 1 and y = L boundaries. The parameter 7 is 
introduced to adjust the variance a v of the velocity field since the (discrete 
time) model is unstable for large values of the variances [0 . 

In our numerical simulations, for all the models we explored, the mean 
scalar profile is linear. If we expand 9 around the mean profile we obtain both 
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an additive noise term, /3v n (i), and a convective (multiplicative) noise term. 
This separation also occurs in the continuum equation. We find that differ- 
ent regimes of behavior of the scalar PDFs can be characterized empirically 
by a single parameter B that measures the relative strengths of convective 
and additive noise terms: B = (To/ where o~q is the variance and the 
characteristic length scale of 9. 

We first consider the regime B < 1 which occurs for a^r c /z/ < 1 where 
r c is the correlation time of the velocity field. Different models are defined 
by the dynamics of V n (i). When the velocity field has Gaussian amplitude 
fluctuations on a time scale r a comparable to the diffusive time scale r d = 1/v 
we obtain PDFs with exponential tails for 9. 

Model A: The easiest way to obtain amplitude fluctuations is to use a 
stochastic model and choose the stream function to be a product of two noise 
terms r} n (i)r}' n (i) where rf n {i) is a white noise and r) n (i) a Gaussian noise with 
a correlation time r a , both having delta-function spatial correlations. The 
stream function and the velocity display two correlation times: (i) the mixing 
time r c governed by r)' n (i) chosen to be one timestep and (ii) the time scale 
r a on which the amplitude varies. This model exhibits exponential tails in 
the PDF of 9(i) when r a r d (See Fig.l). On the other hand when r a « r d 
the distributions are Gaussian. 

Model B: We next provide numerical evidence to show that exponential 
tails also occur when the stream function is derived from a deterministic 
chaotic model provided it leads to amplitude fluctuations on the appropriate 
time scale. We investigate the following model: define the stream function 
in terms of an auxiliary variable <f> n (i) by ip n (i) = V 2 L <p2n(i)', let <f> n (i) be 
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updated according to (f> m +i(i) = \ J2k F((j> m (k)), where the sum includes the 
site i and its nearest neighbours. The function F is chosen to be F(x) = 
(1 — 2a)x+2ax 3 ; we consider a = 2 for which the maximum liapunov exponent 
is 0.95 and the correlation length is measured to be one lattice spacing. 
The PDFs for u,v display Gaussian tails (See Fig.2a). The PDF for 0(i) 
is shown in Fig. 2b which clearly shows exponential tails over five decades 
(similar results persist for a range of values of a). The slope of the tails 
is proportional to y/v/{f3a v ) as expected from dimensional considerations. 
To demonstrate that the velocity field derived from this deterministic model 
displays intermittent amplitude fluctuations we compute the autocorrelation 
functions C 2 (n) =< v n v > — < v > 2 and C 4 (n) =< v^Vq > — < Vq > 2 , 
where < .. > indicates time average. In Fig. 3 we show 64(71) which shows two 
exponential decays associated with time scales r c /2 « 1, r a /2 ps 4. The inset 
shows C 2 (n) which is expected to display a decay time (l/r c + l/r a ) _1 ps 1.6, 
consistent with the above interpretation. Note that in Fig.2b is comparable 
to r a and the kurtosis decreases with increasing r d . 

To elucidate this mechanism we present a mean-field-like toy model with 
only additive noise since the multiplicative noise term is small for B < 1. 
The model is described by the Langevin equations, 

x = -vx + yr}!, (3) 

y = -ay + T] 2 , (4) 

where the variable x couples to a noise with an amplitude y that fluctuates 
in time. The terms rji are assumed to be white noise with < rji >= 
and < r)i(t)r)j(t') >= a?5ij5(t — t'). The crucial time scales are given by 
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r a = 1/at and T4 = 1/v. It is convenient to introduce scaled variables xq = 
xj \fDiD 2l and yo = yj\[D 2 where D\ = o\j2v, D 2 = a 2 /2a. Since yo is 
independent of xq one can obtain a stationary solution for the Fokker-Planck 



(FP) equation |TT| of the form P(x 0l y ) = Q(x \yo)R(yo), where R(yo) = 
[27r] _1 / 2 ea;p(— |/q/2). The conditional probability Q(xo|2/o) satisfies another 
FP equation, that can be solved iteratively for small ajv ||12|| . For large xq 
we have obtained the first three terms of p(xo) = J P(xo,yo)dyo expanded as 
a formal power series: p(x ) = Z)^L (a/z/) n p n (x ). The functions Po,Pi,p 2 
are combinations of modified Bessel functions that have exponential tails for 
large values of xo with pot^o) = [7r]~ ^^i^od^ol), where K is the modified 
Bessel function of order zero. Note that the exponential tail results from 
a linear FP equation due to a branch cut in the characteristic functions in 
contrast to Ref.5 in which a nonlinear FP equation is used with simple poles 
in the characteristic function. In the limit of very large a/u, p(xq) becomes 
Gaussian. We compute the low-order moments and find that the kurtosis is 
given by K = 3 + this interpolates between 9 for small ajv and 3 for 
large ajv corresponding to p and the Gaussian distribution respectively. In 
Fig. 4 we show the PDF of xq obtained from numerical simulations for a = v 
for which K = 6 corresponding to that of an exponential distribution. 

The form of the PDF depends not only on the parameter B but also 
on the nature of the velocity correlations. In addition to exponential tails 
discussed above other non-Gaussian behavior can occur in different models. 
For example, the introduction of another variable z using z = —fiz + xr] 3 
in the toy model (Eq.3,4) leads to a stretched exponential PDF for z with 
an exponent 2/3 for appropriate a, v and fi. A similar generalization of the 
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(stochastic) model A can be constructed. We emphasize that the PDF of the 
variable depends on the distribution of the slowly varying amplitude of the 
noise but not on the distribution of the noise itself. 

We now consider the model in which the velocity is obtained from ip n {i) = 
Vn{i), where r] n (i) is a Gaussian noise correlated over a time r c with a corre- 
lation length £„. The PDFs of 9 are Gaussian when B is smaller than unity 
even for r c pa t&. This is true for both the discrete time and the continuum 
model. However, for B >pa 1 the PDFs are non-Gaussian and we study 
this case using the continuum model i.e., Eq.(l). The resultant PDFs are 
shown in Fig. 5 for different parameters, (i) For B » 1 the distributions 
(Fig.5a,5b,5c) clearly do not have exponential tails; however, strikingly, we 
could fit the PDFs very well with a modified Lorentzian with two parameters 
(k, 5): p(x) — C/(l + kx 2 ) 1+s where C is fixed by normalization. This func- 
tion is a variant of the form proposed by Sinai and Yakhot||. (ii) For B pa 1 
the scalar PDFs appear to have exponential tails over a narrow range of pa- 
rameters (this case corresponds to Fig.5d). However, the fit with a modified 
Lorentzian agrees over a wide range of B values including this narrow range. 

Since the convective noise term is not negligible we model the B >pa 1 
regime by a single variable with both an additive and a multiplicative 
(£2) noise term: 

x = -vx + Ci + (2^- (5) 

When Q are white noise with < Q(t) >= and < d(t)Q(t') >= 2DiSij5(t — 
t'), a straightforward solution of the Fokker-Planck equation yields exactly 
the modified Lorentzian form for the PDF with n = D 2 /Di, 5 = u/2D 2 
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HJ . This provides another example of a simple toy problem that mimics the 
PDFs of passive scalars. 

We note that the exponential tail regime for passive scalars in several 
experiments |2|, [| appears to correspond to the B < 1 regime in our model 
where intermittent velocity fluctuations lead to exponential tails. Experimen- 
tally such amplitude variations may be connected with coherent structures 
such as plumes and turbulent bursts [P^|. A quite different extended dynam- 



ical system in which our mechanism is operative is a Capacitive Josephson 
Junction Array driven by external dc and ac currents; numerical simulations 
in chaotic states indicate exponential tails in the PDFs of local junction 
voltages transverse to the external current |15| . This can be ascribed to in- 



termittent current fluctuations caused by random vortex motion. However, 
establishing the occurrence of such amplitude fluctuations starting from the 
underlying equations, e.g., Navier-Stokes equations with appropriate bound- 
ary conditions, remains a challenge. 
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1 Figure captions 

Fig.l: PDF for the normalized fluctuation X = 59 /ag for Model A {stochas- 
tic) with L = 48, z/ = 0.1,7 = 0.4, (5 = 0.1, r a « 10. The white noise 
terms are uniformly distributed between -1 and 1. The variances are ag m 
0.034, a v m 0.21 and the kurtosis is 4.22. 

Fig. 2: Data for Model B (deterministic). 
a) PDF for v normalized by a v has Gaussian tails, b) PDFs for the normalized 
fluctuation X are shown for two parameters. The dashed line is drawn with 
the same slope as that of the tail, i) Upper curve (shifted up by two decades) 
is for v = 0.1,7 = 0.5, f3 = 0.1. The variances are: ag xs 0.045, a v xs 0.2 and 
the kurtosis is 4.62. ii) Lower curve is for v = 0.05,7 — 0.25, /5 = 0.1. The 
variances are ag ~ 0.031, a v ~ 0.1 and the kurtosis is 4.12. Data are obtained 
using 5 million points in intervals of 10 time steps at a site in the middle of 
a 48 x 48 lattice. 

Fig. 3: Correlation functions of the velocity used in Fig. 2. A semilog plot 
of C 4 (n) vs n is shown. The straight lines have slopes 1.3 and ps 0.24. The 
inset shows a semilog plot of 62(71) vs n; the straight line has a slope 0.69. 
See text for discussion. (The variance a v is set to unity.) 

Fig. 4: PDF for x = x/a x in the toy model of Equations (3) and (4) with 
a — v. The noise variables 771,772 are uniformly distributed with <j\ — 02 — 1 
and the measured a x = 0.5. 

Fig. 5: PDFs for the passive scalar fluctuations with B >k, 1 using the 
continuum model Equation (1). The data are obtained with L = 48, v = 
0.1, P = 0.1 using the discretizations St = 0.05, 51 = 0.5 and velocity correla- 
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tions ~ 51 and r c = 2. The upper three curves (a,b,c) are shifted up from 
the lowest curve (d) by 6,4 and 2 decades respectively. The curves a,b,c,d 
correspond to a v = I.47 where 7 = (0.4,0.35,0.3,0.2) with kurtosis values 
(15.6,6.3,4.1,3.5). The dashed lines are modified Lorentzian (ML) fits for 
the data. 
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